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Abstract— The impulse and square-wave Input response of 
different GaAs metal-semiconductor-metal photodetector (MSM) 
designs are theoretically examined using a two dimensional drift- 
diffusion numerical calculation with a thermionic-field emission 
boundary condition model for the heterojunctions. The rise time 
and the fall time of the output signal current are calculated 
for a simple GaAs, epitaxially grown, MSM device as well as 
for various double-heterostructure barrier devices. The double 
heterostructure devices consist of an AlGaAs layer sandwiched 
between the top GaAs active, absorption layer and the bottom 
GaAs substrate. The effect of the depth of the AlGaAs layer on 
the speed and responsivity of the MSM devices is examined. It 
Is found that there is an optimal depth, at fixed applied bias, of 
the AlGaAs layer within the structure that provides maximum 
responsivity at minimal compromise in speed. 


I. Introduction 

M ETAL-scmiconductor-mctal (MSM) photodetectors are 
becoming increasingly attractive in optoelectronic com- 
munication systems, high-speed chip-to-chip connections, and 
high-speed sampling applications [1]. The principal advantages 
of MSM photodetectors that make them an excellent choice 
for on-chip detectors are their responsivity-bandwidth perfor- 
mance, compatibility with existing planar integrated circuit 
technologies, and relatively low voltage operation [1J. Im- 
provement in MSM detector design structures, specifically the 
introduction of a GaAs/AlGaAs heterostructure, has advanced 
the realization of these applications [2]. 

The basic MSM detector structure is shown in Fig. 1. This 
structure consists of interdigitated metal fingers formed on the 
top surface of a semiconductor layer. Light incident on the top 
surface of the MSM structure is absorbed within the underlying 
semiconductor resulting in the creation of electron-hole-pairs 
(EHP's). The application of a bias to the metallic fingers 
creates an electric field within the underlying semiconductor 
which acts to sweep the photogenerated carriers out of the 
device. How fast these carriers are collected and how many 
of them actually survive to the contacts within a particular 
collection time determine the speed and the responsivity, 
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Fig. 1. Schematic diagram of an interdigitated metal-semiconductor-metal 
(MSM) structure. 

respectively, of the detector. Carriers generated deep within the 
semiconductor must traverse a greater distance before they are 
collected at the contacts compared to those generated near the 
surface. Depending upon the magnitude of the electric field 
within the semiconductor, the time needed to collect those 
earners generated deep in the device can vary drastically. 
Under low bias conditions, ^magnitude of 5-10 V, which is 
typical for most integrated circuit applications, this collection 
time can be prohibitively long in high speed applications. 

It has been suggested that the insertion of a double het- 
erostructure layer can improve the time response of MSM 
detectors 12]. The double heterostructure layer acts to block 
those carriers generated deep within the device structure. As 
a result, only those carriers photogenerated within the top 
absorption layer are collected leading to a fast overall response. 
However, as described above, the responsivity principally 
depends upon the number of photogenerated carriers collected 
at the contacts. A high responsivity, especially at low input 
power levels, dictates that most of the photogenerated earners 
be collected. Since many of the photogenerated carriers are 
produced deep within the semiconductor layer, the insertion of 
a double heterostructure layer in order to improve the speed 
of the device, necessarily reduces its responsitivity as well. 
Therefore, there exists a fundamental tradeoff between the 
speed of response and responsivity of a heterostructure MSM 
detector. 

It should be noted that the absorption coefficient varies 
strongly with the optical wavelength. In this paper, we ex- 
amine the response of the photodetector to only one incident 
wavelength, 840 nm. The speed of response of the detector can 


00 1 8-9383/94S04.00 C 1994 IEEE 




SALEM et at.: THEORETICAL STUDY OF THE EFFECT OF PHOTODETECTOR PERFORMANCE 


1113 


be quite different at other wavelengths since the absorption 
may occur closer to the surface or deeper into the device. 
If for example, the absorption depth is much smaller than 
the interdigitated finger spacing then the detector’s speed of 
response is limited by the finger separation and the placement 
of the AlGaAs layer has minimal effect. However, at an 
incident wavelength of 840 nm the absorption depth and finger 
separation are comparable and the placement of the AlGaAs 
layer becomes important. Therefore, the tradeoff between the 
speed of response and the rcsponsitivity discussed herein 
applies principally to incident light at 840 nm. 

In this paper, we present results of a theoretical study of 
the performance of a GaAs based MSM detector examining 
the effect of including a sandwiched AlGaAs layer between 
the active GaAs layer and the substrate. A two-dimensional 
drift-diffusion numerical calculation in combination with a 
thermionic-field-cmission model is employed in this study. 
The drift-diffusion/thermionic emission model along with the 
details of the numerical solution method used are described in 
Section D. Calibration of the model and comparison to existing 
models is discussed in Section III. In Section IV different 
GaAs MSM double heterostructure devices are examined using 
the model. In Section IV we examine the effect of the AlGaAs 
layer on the rise time, fall time, and the magnitude of the 
output signal current and optimize the placement of such a 
layer as it relates to these parameters. 

n. Model Description 

To understand the behavior of photo-generated carriers 
under the influence of drift and diffusion forces, Poisson’s 
equation, the current continuity equations for electron and 
holes, and a rate equation for charged traps are solved on a two 
dimensional simulation domain for the electrostatic potential, 
and the carrier concentrations. The drift -diffusion approach 
is sufficient to characterize the operation of MSM detectors 
since these devices are typically of the order of microns in 
length and width. In this study, the detectors examined are 
restricted in dimensions to several microns where the drift- 
diffusion approximation is valid. The basic equations used in 


the simulator are [3], [4]: 

e 0 £ # VV = — q(p - n + Nd — Na + ***) 0) 

-±VJ n -G n + R n + ^= 0 ( 2 ) 

q Ot 

±VJ p -G p + R t + ^ = 0 ( 3 ) 

J n = + qD n Vn (4) 
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“J7 = (Rn - C n )sRH - (Rp - Gp)sRH (6) 
at 

where is the electrostatic potential, n and p are the electron 
and hole carrier concentrations, Nd,a are the donor/acceptor 
concentrations, fit is the electron-filled trap concentration, 
J np are the electron/hole currents, Rn tP and G njt are the 
electron/hole recombination and generation rates, /i„ iP and 
D n>p are electron/hole mobilities and diffusivities, and the 
subscript SRH in (6) denotes the Shockley-Read-Hall events. 


In most situations, the trap density is such that the effects 
of (6) are negligible [3]. The terms Rn %p used in (2) and 
(3) represent the sum of Shockley-Read-Hall, radiative, and 
Auger recombinations for electrons and holes 14]. The optical 
generation rate, G, is generally assumed to be gaussian for the 
impulse response and is typically given as [5] 
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where is the peak optical power density incident on 
the surface, hv is the photon energy, <**b» is the absorption 
coefficient, x„ is the location of the center of the beam, t 0 is the 
location of the optical input peak, <r t and o x are related to the 
full width, half maximum (FWHM) of the pulse. In this paper, 
we examine the response of only a representative unit cell of 
the device in which the illumination is practically uniform. ' 
The beam is assumed to be centered within the unit cell. Given 
the dimensions of the unit cell, the lateral decay of the beam 
is negligible within the calculations presented here. Therefore, 
for simplicity, uniform illumination is assumed laterally within 
the unit cell and the spatial exponential term in (7) is set 
to unity. The metallic fingers are assumed to be completely 
transparent In practice, there are of course shadowing and 
reflective losses at the metallic surface which would lead to a 
reduction in the amount of photogenerated carriers within the 
underlying semiconductor material. For simplicity and due to 
a lack of detailed information about the extent of these losses, 
we assume here that no losses occur. 

A standard field-dependent mobility for the electron mobil- 
ity, /A. 14] is used. 


Mn = 


F 3 

Mno "1" pi 



( 8 ) 


where Uno is the zero field mobility, F a is the critical electric 
field, F is the local electric field, and v, is the saturation veloc- 
ity. The hole mobility is assumed to be constant, independent 
of the electric field. It is further assumed that the Einstein 
relation holds for the diffusivities. The use of the Einstein 
relation, though universally accepted, is not totally correct 
since it applies only strictly to equilibrium. Improvement over 
the Einstein relation necessitates determining the mobility 
and diffusivity using more exact methods of solution of the 
Boltzmann equation, i.e., the ensemble Monte Carlo technique. 
Presently, this information is not available to us and thus for 
simplicity, the Einstein relation is adopted as is typically done 
in drift-diffusion solutions. 

The model’s modifications made to handle heterojunctions 
follow the same approach as given by Sutherland and Hauser 
(6]. In addition, a thermionic-field emission boundary condi- 
tion is used to specify the current density at heterojunction 
interfaces. The thermionic emission boundary condition [7], 
{8] is implemented in parallel with the drift-diffusion model. 
The actual current across the heterointerface is limited by 
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either the thermionic emission or diffusion current, depending 
upon which is smaller. In the present model, die drift-diffusion 
and thermionic emission current densities at the heterointerface 
aie computed and compared. The actual current density across 
the heterointerface is then given by the smaller of these two 
quantities ensuring that the current across the junction is either 
thermionic emission or diffusion limited. Tunneling across the 
heterojunction is modeled based on the theory of Crowell and 
Rideout [9] and as developed in [10]. The tunneling current 
density is added at the heterojunction interface to give the total 
current density across the interface. 

Equations (1 )-(6) along with the proper boundary conditions 
can be solved for the three fundamental variables S', n, and p. 
The Dirichlet boundary conditions of the electrostatic potential 
used at the Schottky contacts are 

= fa + p P - il>» (9) 

where is the built-in potential, * app is the applied bias 
voltage, and is the Schottky barrier height. From the 
thermionic emission and diffusion theory of Crowell and Sze 
[11] the carrier concentrations at the Schottky contacts are 
specified in terms of the current density passing through them 
as [3], [4], [12]: 


Jn * n = -qv n (n - n 0 ) 

J T • h = qv p (p - p 0 ) ( 10 ) 

where t>„ iP are the electron/hole thermionic recombination 
velocities, n is the unit normal vector, and n 0 and p 0 are 
the equilibrium electron and hole carrier recombination at the 
Schottky contacts. Specifically, n Q and p Q are given as 

n 0 =N c e ( -^' ksT > ( 11 ) 

p Q = fq ve (9(^-E^p)/kBT) ( 12 ) 


where N c and N v are the electron and hole effective density of 
states, and is the energy gap. At any interface, Gauss's 
law can be applied to relate the normal component of the 
electric flux density to the interface charge as 


dyp 
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where Qi n t is the interface charge. The free surface carrier 
concentrations are determined from conditions on the normal 
current [4], [12]: 


J n -n=q - i?* urf (14) 

J p n=ql?'" f (15) 

where iJ* urf is the surface recombination rate. Equations 
(12>-{14) reduce to the usual Neumann boundary conditions 
when the interface charge and the surface recombination rates 
vanish. 

The box integration approach of the finite difference method 
[4] is used to discretize ( 1 )— <6). These equations are then 
solved with their appropriate boundary conditions on a nonuni- 
form, two-dimensional mesh consisting generally of 75 x 33 
points. The Scharfetter-Gummel technique [4], [13] is applied 
in the formulation of the discretization equations. The resulting 


system of equations is linearized using Newton's method. 
Due to storage and computation time limitations the Gaussian 
elimination method is deemed inappropriate for the solution 
of the linearized system. Therefore, the system of equations is 
solved using an iterative approach known as the bi -conjugate 
gradient squared (BICGS) method [14], [15]. Only nonzero 
elements in the matrix are stored using this approach greatly 
reducing the data storage requirements of die code. After 
the program reaches convergence, the potential and carrier 
concentration profiles are readily available, from which most 
of the macroscopic variables of interest can be calculated. The 
current densities, J n and J p , are calculated from (4) and (5) 
using Scharfetter-Gummel ’s exponential scheme [13] for the 
carrier concentrations. 


ID. Calibration of toe Model 

To demonstrate the accuracy of the model described above 
it may seem at first that direct comparison to experimental 
measurements would be best. However, this is not an easy 
task since information about the experimental setup is actu- 
ally needed to accurately compare the theoretical results to 
experimental measurements. Different authors [5], [16] have 
attempted comparing the calculated response of GaAs MSM 
photodetectors from the drift diffusion model to experimental 
measurements. Landheer et al [16] found that an equivalent 
circuit model for the experimental setup is needed to accurately 
compare the theoretical results to experimental measurements. 
The output current predicted by the circuit model of Landheer 
et al [16] did not lead to very accurate agreement with the 
experimental measurements. This discrepancy is apparently 
due to the simplistic treatment of the external circuit and the 
lack of complete knowledge of the experimental circuit param- 
eters. Alternatively, Sano [5] proposed an analytical model for 
GaAs MSM photodetectors based on the solution of the drift 
diffusion equations (1 )— (6). The analytical model used is an 
equivalent RC circuit model in which the resistance, R , and 
the capacitance, C, are calculated based on the electric field 
and carrier concentrations obtained from the drift diffusion 
solution. The equivalent circuit model was implemented in a 
SPICE-like circuit simulator and the transient responses from 
the circuit simulation were compared with measured responses 
for an MSM photodetector. Though Sano [5] achieved good 
agreement between the circuit simulator and the experimental 
measurements, this is accomplished by numerically adjusting 
die circuit parameters in the model to ensure agreement with 
the experiment. Therefore, direct comparison of the drift 
diffusion results to experiment is hampered by the insufficient 
information about the experimental setup used in the measure- 
ments. In the absence of information about the experimental 
setup we choose to compare our model instead to other existing 
numerical models. 

We have compared the results obtained using the present 
model against several existing device simulators and against a 
1-D analytical solution for a simple 1 /im GaAs p-n junction 
diode. The donor and acceptor concentrations chosen for this 
test device are Nd = N& = 10 17 cm“ 3 and the carrier 
mobilities are /in = 7000cm 2 /V-s and - 300cm 2 /V-s. 
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Fig. 2. Comparison of the electrostatic potentials obtained from three differ- 
ent numerical solutions for a GaAs pn junction diode forward biased at 1 V: 
cur r e nt program (circles), PC- ID (diamonds), and STEBS-2D (squares). 



Fig. 3. Comparison of the electron and hole concentrations for a GaAs pn 
junction diode forward biased at 1 V calculated from the current program 
(circles), PC-1D (diamonds), and STEBS-2D (squares). 

Fig. 2 shows the calculated electrostatic potential V obtained 
from the present model and from two other simulators, PC- ID 
[17] and STEBS-2D [18]. PC- 1 D is a commercial, one- 
dimensional, drift-diffusion model released by Sandia National 
Laboratory and STEBS-2D is a full hydrodynamic model 
developed at the Georgia Institute of Technology. The car- 
rier concentrations obtained from the three models are also 
presented in Fig. 3. As can be seen from inspection of Figs. 2 
and 3 , the present program shows precise agreement with the 
other two models under comparable conditions. 

The electrostatic potential within the p-n junction diode 
calculated using the present model can also be compared 
to a simple one-dimensional analytical solution of Poisson’s 
equation using the depletion approximation. Inspection of 
Fig. 4 shows good agreement between the numerical model 
and the analytical model except near the edge of the depletion 
region where the depletion approximation fails to properly 
account for the tails in the carrier distributions. 

IV. Application to MSM Photodetectors 

The present numerical model is applied to study the re- 
sponse of GaAs based metal-semiconductor-metal (MSM) 
photodetectors. The general device structure considered is 
sketched in Fig. 1 above. The barrier height assumed for 
the Schottky contacts in the structure is 0.7 V and the 
GaAs material is taken as semi-insulating with a background 
doping concentration of Np = 10 15 cm~ 3 . The interdigitated 



Fig. 4. Comparison of the electrostatic potentials for a GaAs pn junction 
diode forward-biased at 1 V calculated from the present model (circles) 
and a one ■dimensional analytical solution (diamonds) using the depletion 
approximation. 


TABLE 1 

Material Parameters Used 


Parameter 

Units 

GaAs 

Ref. 

AlGaAs 

Ref 

Electron Mass (m* ) 

— 

0.061 

(18] 

0.088 

(18) 

Hole Mass (mj) 

— 

0.48 

(18] 

036 

(18] 

Dielectric Constant (f ) 

— 

13.1 

(18] 

12.4 

(18) 

Energy Gap (E p ) 

(cV) 

1.42 

118] 

1.74 

(18] 

Electron Affinity (x) 

(eV) 

4.06 

[19] 

3.80 

(19] 

Electron Mobility (pn) (an 2 /Vs) 

7000 

119] 

2500 

(19] 

Hole Mobility (p p ) 

(an* /Vs) 

300 

119] 

150 

[19] 

Electron Lifetime (r„ ) 

(*) 

10 -T 

(20] 

10 -i 

[19) 

Hole Lifetime (r p ) 

(I) 

lO -7 

(20] 

10 -8 

[19] 

Radiative coef. (B) 

(cm»/s) 

2.04 xlO -10 

(21] 2.04X10 -10 

a 

Electron Auger Coef 

(cm 6 /*) 

1.6 x 10 -s ® 

(21) 1.6 x 10 -M 

a 

(Cn) 






Hole Auger Coef. (C p ) (cm e /s) 

4.64 XlO -29 

[21] 4.64 xlCT” 

a 

Absorption Coef. 

(cm -1 ) 

10* 

(22] 

0 

b 

(o.b.) 







■Due to lack of information for AlGaAs, GaAs values arc used. 
b No absorption takes place for the AlGaAs at A as 0.84 /im. 


TABLE D 

Input Signal Parameters 


Peak Power (U p ) 

0 68 mW/cm 2 

Wavelength (A) 

840 nm 

Beam Diameter (D) 

60 pm 

Beam’s Peak Position (A”©) 

23 pm 

Signal Peak (t©) 

10 ps 

Full Width Half Max. (FWHM) 




metallic finger widths and spacings are 1 jim and 3 ixm, 
Tespectively. The material parameters used in the calculation, 
i.e., the zero field carrier mobilities, lifetimes, etc. art compiled 
in Table I. These parameters are compiled from [19H23]. 
The parameters used to characterize the optical input signal 
axe collected in Table II. In these calculations the boundary 
condition at the bottom surface and along the sides of the 
device is assumed to be floating. 

The response of the GaAs MSM device as sketched in 
Fig. 1 is compared to a double heterostructure device. The 
double heterostructure device consists of an AlGaAs barrier 
layer of thickness d 2 , sandwiched between two different GaAs 
layers, of thicknesses di and ^ 3 , as shown in Fig. 5. The 
Al concentration within the AlGaAs layer is assumed to be 
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X-Axis 

Fig. 5. Schematic diagram of the double-barrier heterostructure MSM show- 
ing the sandwiched AlGaAs layer. The solid shaded regions on the lop of the 
figure represent the Schotiky contacts. 

held fixed at 25%. Different double heterostructure devices are 
examined with variable depths of the AlGaAs barrier layer. 
As a result, the active, photo-absorption layer thickness du 
obviously also varies in these devices. 

As described above, the primary function of the buried 
heterostructure layers is to block the collection of the carriers 
photogenerated deep within the GaAs substrate by preventing 
them from diffusing back into the active layer and towards the 
contacts. In addition, at low applied bias, the barrier between 
the active and the AlGaAs layers acts to confine the photogen- 
erated carriers within the active region. At higher applied bias, 
the heterostructure blockage of the photogenerated carriers 
becomes less effective due to the much greater band bending 
present in die device. Subsequently, the location of the AlGaAs 
layer greatly effects the charge collection attributes of the 
device depending upon the field distribution and the applied 
bias. In the structures examined here, the background doping 
is n-type, implying that the primary photogenerated carriers 
collected are holes. 

The calculated impulse response at different active layer 
thicknesses, along with the corresponding rise and fall times of 
the signal are displayed in Fig. 6. Five different structures are 
examined, four double heterostructure devices and one GaAs 
bulk device 6 /im in thickness. The double heterostructure 
devices consist of a top GaAs layer, d\ % ranging in thickness 
from 0.5 to 4 /xm, a 1 /xm AlGaAs layer, d 2 , followed by 
a GaAs epilayer, d$, ranging in thickness from 4.5 to 1 /xm. 
Notice that the total width of all three layers combined remains 
constant at 6 /im. The rise and fall times are defined as the 
time it takes the output signal to go from 10% to 90% and 
from 90% to 10% of its maximum value, respectively. The 
voltage applied to the device is -5 V. As can be seen from 
Fig. 6, the fastest response occurs for the device configuration 
with a 0.5 /xm active layer thickness. This is obvious from 
both the curve corresponding to the 0.5 /xm device as well 
as from its corresponding fall rime. However, the maximum 
output signal magnitude for the 0.5 /xm device is significantly 
less than for the other cases. This is as expected, since the 
response speed is achieved at the expense of lower output 
signal magnitude because the slower earners, those generated 
deep within the device, are blocked from being collected by 
the heterojunction barrier. As the active layer thickness d\ 
increases to 1 /xm, more carriers are generated within the 
top, active GaAs region. As a result, a greater number of 



Fig. 6. The calculated impulse response at different active layer thicknesses 
for the double-barrier MSM heterostructure at an applied bias of -5 V. 
The doping levels are 10 I5 cm“ 3 for the GaAs layers and 10 1 axT 
for the AlGaAs layer. Curve 5 is for the case without the AlGaAs layer. 
The numbers in the legend re p rese nt the layer thickness in microns of the 
GaAs/AJGaAs/GaAs layers as shown in Fig. 5. Also included are the rise and 
fall times for each of the 5 curves. 
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Fig. 7. The maximum output current signal as a function of active layer 
thickness for the double-barrier MSM beterostnicture. 

photogenerated carriers are collected producing a higher output 
signal current Though the d\ = 1 /xm device does not show 
as rapid a collection of the photogenerated carriers as the 
d\ - 0.5 /xm device, the field is sufficiently strong and the 
carriers are still relatively close to the collecting contacts that 
a reasonably high speed of response is retained; a fall rime of 
22 ps is achieved as compared to 13 ps for the 0.5 /xm device. 
If the active layer thickness is increased further to 2 /xm, the 
same trend is observed; more carries are collected from the 
bottom of die active layer resulting in a longer fall rime, ^42 
ps, and slower speed of response. Interestingly, the maximum 
output signal current ultimately decreases with increasing d\. 
The maximum signal current is plotted as a function of active 
layer thickness in Fig. 7. As can be seen from Fig. 7, the output 
signal current reaches a maximum for an active layer thickness 
of d\ = 2/xm. The maximum signal current is significantly 
less for a device with d\ = 4 /xm than with d\ — 2 /xm, though 
the fall times are comparable. Clearly, at an applied bias of 
-5 V t there exists an optimal thickness, in terms of speed of 
response and collection efficiency, for layer d\. 

In an attempt to understand the origin of the peak in the 
output current signal versus active layer thickness curve, we 
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fig. 8. Electric field vector plot at the time-step corresponding to the peak 
current showing the direction of the electric field at each mesh point within the 
two-dimensional simulation grid used for the 2.0 pm active layer thickness 
device. The AJGaAs layer is located 2.0 //m from the lop and is 1.0 /im in 
thickness. Only the top 3 of the device is shown in the figure. The right lop 
contact of the device is biased at —5 V with respect to the left top contact The 
horizontal channeling of the electric field at the middle of the heterobarrier 
is due to the accumulation of electrons within the AJGaAs trapped by the 
potential band bending at either interface. Notice that the lateral component 
of the electric field vanishes within the active region above the heterobairier. 


fig. 9. Electric field vector plot at the time-step corresponding to the peak 
current showing the direction of the electric field at each mesh point within the 
two-dimensional simulation grid used for the 4.0 jim active layer thickness 
device. The AJGaAs layer is located at 4.0 pm from the top and is 1.0 
pm in thickness. Only the top 3 pm of the device is shown in the figure. 
The right top contact of the device is biased ai —5 V with respect to the 
left top contact Notice that the electric field “swirls” around deeper within 
the active region of the device. The presence of the lateral component of 
the electric field acts to increase the trajectory of the carriers leading to a 
longer time between generation and collection with the subsequent increase 
of recombination losses. 


have examined the electric field profiles within each device. 
Field vector plots showing the direction of the electric field 
at each mesh point within the device are shown in Figs. 
8 and 9 for the 2.0 and 4.0 /xm active layer thickness 
devices, respectively. Comparison of the two figures clearly 
shows that the lateral field component essentially vanishes 
near the heterostmeture in the 2.0 /xm width device, while 
a significant lateral component persists within the 4.0 /xm 
width device. The electric field points vertically away from 
the heterostructure acting to accelerate the photogenerated 
holes towards the collecting contacts in the 2.0 /xm device. 
Alternatively, in the 4.0 pm device, the lateral component 
of the electric field acts to acclerate the holes to some extent 
laterally. As a result, the hole trajectory for collection is longer, 
effectively delaying their collection at the contacts. Since the 
holes spend more time in this region their chances of suffering 
a recombination event increase, leading to a reduction in 
the collected current. Subsequently, this could result in an 
increased collected current within the 2.0 /xm width device 
than in the 4.0 ixm width structure. A more detailed analysis of 
this problem using a more sophisticated, hydrodynamic model 
will be made in the future to further test this hypothesis. 

The effect of the heterojunction on the response to a square 
wave input is considered next. The same model and approach 
as described above is used but with a square wave input 
signal applied to the device. The output current response to 
a square wave input with repetition rate of 1 GHz and 0.5 ns 
duration time is shown in Fig. 10 for different GaAs active 


layer thicknesses. The applied voltage is again -5 V. It is 
clear from this figure that the output signal amplitude peaks 
at an active layer thickness of 2 pm and then decreases as the 
active layer gets thicker. The explanation for this is similar to 
that described above for the impulse response. The rise and 
fall times however are not quite consistent with the impulse 
response results. This is due primarily to the definition of these 
quantities. Although the 0.5 /xm curve in Fig. 10 seems to 
decay faster than any of the other cases, the fall time, for 
example, is not the shortest. This is due to the fact that the 
90-10% measure of the fall time depends on the magnitude 
of the signal peak itself. Since the signal peak is very much 
•mailer in the 0.5 /xm device, its decay to 10% takes longer 
than that for a much higher signal peak, though the signal in 
the 0.5 /xm device has decayed quicker below some threshold 
level. Subsequently, the definition of the rise and fall times 
typically given, 10-90% and 90-10% respectively, are of 
questionable value when evaluating the performance of an 
MSM detector subject to a square wave input. However, it is 
clear that the presence and the location of the AlGaAs layer is 
critical to the performance of these photodetectors and that the 
optimal location of the heterojunction barriers is a function of 
the applied bias. This is clearly demonstrated in Fig. 1 1 where 
the applied voltage is increased to -10 V. We note that in this 
case the peak in the collected current shifts to the 3 /xm active 
layer thickness device. The shift in the peak collected current 
occurs since the field depletes deeper into the semiconductor 
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Fig. 10. The calculated output current response of the double-bairier MSM 
beterostructure to a square-wave input at different active Layer thicknesses and 
under an applied voltage of —5 V. Curve 6 it for the case without the AIGaAs 
layer. The numbers in the legend rep rese nt the Layer thickness in microns for 
the GaAs/AIGaAs/GaAs layers as shown in Fig. 5. Alio included are the rise 
and fall times in picoseconds for each of the curves. 



Fig. 11. The calculated output current response of the double-bairier MSM 
heterostructure to a square- wave input at different active layer thicknesses and 
under an applied bias of —10 V. Curve 6 is for the case without the AIGaAs 
layer. The numbers in the legend represent the layer thickness in microns of 
the GaAs/AIGaAs/GaAs Layers as shown in Fig. 5. Also included are the rise 
and fall times in picoseconds for each of the curves. 

material. Subsequently, more carriers are swept out by the 
relatively high-field producing a greater collected cunent at 
die contacts. 

Due to computational limitations, it is presently possible 
to simulate the entire substrate. Subsequently, only a small 
portion of the substrate can be simulated. In this case, the 
bottom surface can be treated either as floating or as an ohmic 
contact. The use of an ohmic boundary condition along the 
bottom surface leads to qualitatively similar results as those 
with the floating surface condition but with some quantitative 
difference. The appropriateness of these two boundary condi- 
tions needs to be further investigated, which wil be reported 
in a future work. 




Fig. 12. TWo-dtroensiooal conduction band diagram of the double -barrier 
MSM beteroatructure for the 2 pm active layer thickness device at an applied 
bias of -5 V. 

Finally, we examine the dark current for die device sketched 
in Fig. 1 assuming a finger length of 100 pm, typical for most 
MSM structures. The dark cunent is found to be on the order 
of 0.6 pA. Moreover, it is also found that it is practically 
insensitive to die applied bias in the range from zero to —20 
V. This is because, at a 3 pm finger spacing, and at the 
applied biases considered here, the built-in potential barrier 
at the grounded Schottky contact persists as shown in the 
two-dimensional conduction band plot for the 2.0 pm active 
layer thickness device in Fig. 12. However, if die voltage 
is increased or if the finger spacing is reduced, the applied 
voltage at the biased contact punches through to the grounded 
contact reducing the built-in potential barrier. As a result, an 
increased electron dark cunent develops. It is expected then 
that for a smaller finger spacing (~1 pm) and a doping level of 
10 16 cm' 3 that die dark current will increase with increasing 
applied voltage [24], [25]. In that case, a top AIGaAs layer 
[26] can play an important role in limiting the dark cunent 
of the device. 

V. Conclusions 

We have theoretically investigated the effect of including 
a double heterostmcture barrier on the response of MSM 
interdigitated photodetectors using a two-dimensional drift- 
diffusion model with a thermionic-field emission boundary 
condition. To establish the validity of the model we have 
compared it to existing models and to a one-dimensional 
analytical solution. Excellent agreement with these models 
is obtained. The drift-diffusion model is then used to study 
the effect of the presence of a sandwiched AIGaAs layer on 
the response, as measured by the responsivity and speed, of 
MSM photodetectors. The AIGaAs layer introduces a double 
heterostmcture barrier whose location in the structure greatly 
affects the movement of the carriers. We have found that there 
exists an optimal location for these heterojunction barriers 
for which both high speed and high responsivity can be 
achieved for an input optical signal of 840 nm wavelength. The 
placement of the AIGaAs barrier depends on the applications 
in which the MSM’s are to be used. In high speed applications, 
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and under low bias conditions, the AlGaAs layer should be 
placed near the top surface of the device. To obtain a fast 
response under low light illumination levels, the active layer 
thickness and the applied bias should be increased. However, 
operation at high bias can increase the dark current levels of 
tire device if the electrode spacing is small. Subsequently, there 
exist several tradeoffs in the design of an interdigitated MSM 
photodetector. Use of a simulator such as the one described 
here, is essential for optimizing a structure for a particular 
application. 
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